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ABSTRACT 



The radiative and spatial characteristics of 63 ship tracks are analyzed in cloudiness transition 
regions using AVHRR satellite data. Channels 1, 2, 3, 4 and 5 are utilized to determine the 
variations in ship track albebo and temperature in three distinct cloud environments: stratus, broken 
stratus and cumulus. Stratus Tracks in a stratus environment are brightest with respect to broken 
stratus and cumulus tracks. Cumulus tracks have the highest percentage increase in albedo over their 
environment in Channel 1, while broken stratus tracks have the highest percentage change over their 
envirorunent in Channel 3. There is little temperature difference between the ambient cloud and ship 
track in all three environments. Stratus track tend to be longer than either broken stratus or cumulus 
tracks. No significant difference in the width of the track was observed between the environments. 
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I. 



INTRODUCTION 



A. OVERVIEW 

Ship tracks have been observed in satellite imagery for 
many years. In the mid 1960's Conover (1966) described 
"anomalous cloud lines" and attributed their formation to 
passing ships. However, it is only in recent years that 
research efforts have been directed towards explaining the 
atmospheric conditions responsible and cloud physics required 
for ship track formation. These research results have many 
applications ranging from the effect of anthropogenic aerosols 
clouds, which contributes to our understanding of the global 
climate, to intelligence applications, which focus on tracking 
and classifying ships based their signature in the ambient 
cloud field. Ship tracks form in a wide variety of marine 
stratiform clouds, transitioning from uniform solid stratus, 
broken stratus to open- celled cumulus. Each of these cloud 
types have unique formation processes, micro-physical 
properties and radiative characteristics. Since ship tracks 
have not been observed in completely clear air regions it is 
apparent that a link exists between the formation conditions 
necessary for stratiform clouds and those of ship tracks. 
This link will be investigated by describing the spatial and 
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radiative characteristics of the ship tracks which form in 
marine stratiform clouds, specifically the cloud forms 
associated with transitional regions between broken or partial 
cloud and solid stratiform cloud decks. The purpose of this 
study is to determine if the observed variations 'in the cloud 
forms found in these transition regions are reflected in the 
ship tracks which form in them. 

B . BACKGROUND 

1. Ship Tracks 

Ship tracks are formed primarily during the sxammer 
months in the low stratus and stratocumulus clouds of eastern 
ocean basins. Current theory suggests they are the direct 
result of a ship passing under a stratiform cloud and either 
modifying the existing cloud or creating new cloud by 
expelling exhaust gases into the saturated marine boundary 
layer. These gases contain heat, moisture and combustion 
particles, the last of which serve as cloud condensation 
nuclei (CCN) for water droplets. The injection of additional 
CCN into the stratiform cloud deck, by direct means or by gas- 
to-particle conversion of the exhaust gases (Hobbs, et 
al,1980), creates cloud droplets with much smaller radii than 
those found in the adjacent cloud. Since the concentration and 
size of cloud droplets is determined primarily by the 
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concentration of CCN in the air (Radke, et al 1989) , any 
increase in the number of CCN or decrease in cloud drop radii, 
changes the micro-physical characteristics of the cloud. This 
shift in cloud drop properties changes the reflectivity of the 
modified cloud making it readily discernable at certain 
wavelengths in satellite imagery. A comparison of the radiance 
in the NOAA-10 Advanced Very High Resolution Radiometer 
(AVHRR) near infrared channel 3, centered at 3.7 /xm, with the 
visible channel 1, centered at .63 /xm, clearly reveals these 
features. Figure 1. 

Additionally, the buoyancy flux generated by the 
temperature difference between the exhaust gases and 
atmosphere enhances mixing. This buoyancy flux induced mixing 
in concert with the small-scale vortex created by the ship 
motion may initiate and further amplify the ship track cloud 
formation process by introducing marine aerosols from the 
surface of the ocean into the well mixed marine boundary layer 
(Figure 2) . Continuing gas- to-particle conversion of the 
exhaust pltime coupled with the tendency for smaller droplets 
to remain in the cloud longer helps explain the longevity of 
these tracks. A more detailed discussion of proposed formation 
mechanisms and the radiative signature of ship tracks is 
provided by Hindman, 1990 and Coakley, et al., 1987, respectively. 
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Figure 1. Ship Tracks in Satellite Imagery at .63 /xm, (top), 
and 3 . 7 (bottom) . 
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Ship-passage Induced Mixing. 



2 . Marine Stratiform Clouds 

Marine stratiform clouds and their associated 
environments have been the focus of numerous research efforts 
in recent years because of their impact on the planet's albedo 
and global climate. The physical mechanisms which control the 
formation of these clouds are generally explained, but the in- 
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cloud mechanisms responsible for their modification, breakup 
and dissipation are not as well understood. 



Certain 


atmospheric conditions 


must be 


present 


for 


stratiform 


clouds to develop. In 


general , 


stratus 


and 



stratocTomulus clouds form in the summer in the eastern oceans 
when the subtropical high is fully developed. Strong 
subsidence, associated with the subtropical high creates a 
strong temperature inversion over the entire region which 
keeps the moisture confined near the surface (Brest et al, 
1982) . Northerly winds along the coast cause the upwelling of 
deeper, colder ocean water that cools the marine boundary 
layer to saturation. This provides an ideal environment for 
the formation of statiform clouds. 

Synoptically , this inversion slopes from east to west with 
the strong, low inversion in the east near the coast and the 
weak, higher inversion to the west. This occurs on a scale 
comparable to the region between California and Hawaii in the 
eastern Pacific Ocean. Associated with this sloping and 
weakening inversion is a corresponding transition in the cloud 
types from low stratiform clouds in the east to cumuliform 
clouds in the west. 
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On the mesoscale, the reverse has been observed. Betts and 
Boers (1990) and Paluch and Lenschow (1991) observed that 
cloudiness transition regions exist in the marine boundary- 
layer off the west coast of southern California. They are 
characterized by a gradual change from a clear boundary layer 
through small cumulus and broken stratocumulus to a deck of 
solid stratocumulus (Betts and Boers, 1990) . Within these 
cloudiness transition regions the inversion height changes 
from weaker and somewhat lower (950 mb) in the cloud free 
region to relatively stronger and higher (935 mb) in the solid 
stratus clouds (Figure 3) . 

Both Cahalan and Snider (1989) and Paluch and Lenschow 
(1991) noted increases in the liquid water content (LWC) , 
mixing ratio and vertical motion moving from clear air to 
solid stratiform clouds in similar transition regions. These 
variations manifest themselves as changes in the radiative 
signatures of the clouds. 

The structure of the clouds depends largely on the 
environment in which they form. Uniform stratus is 
characterized by a fairly flat top surface which reflects the 
presence of a very strong inversion which suppresses vertical 
motion at the cloud top. A fairly new cloud will have a nearly 
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Figure 3. Mean Profiles for Three Regimes: Cumulus, Broken, 
and Stratus Clouds (Betts and Boers, 1990) . 



homogenous internal structure having not been s\ibjected to 
differential heat fluxes from both the surface and from above. 
As the cloud gets older and thicker through mixing, 
entrainment and a variety of heat fluxes, its texture changes. 
Differential heating allows vertical plumes to penetrate the 
inversion and distort the top of the cloud. As a result it is 
transformed into a solid field of symmetric cellular clouds. 
In the presence of surface winds speeds greater than 10 ms'^ 
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this same cloud field can be modified into cloud streets which 
align themselves parallel to the wind direction. These cloud 
streets consist of the same cellular cloud field but they are 
elongated by the wind. A secondary roll circulation can form 
within the cloud. In between the cloud streets subsidence from 
the secondary circulation creates a variegated cloud top 
appearance. Further entrainment and heat fluxes can cause the 
solid cloud deck to break up producing a cloud field with a 
broken or patchy appearance. If the cloud layer reaches the 
precipitation stage, evaporation of the precipitation moistens 
and cools the air below. Eventually this causes cumuli to form 
below the stratus deck. In time the stratus clouds can 
dissipate due to further entrainment leaving behind only the 
cumulus clouds (Paluch and Lenschow, 1991) . This cloud type 
can be modified by the ship or provide a suitable environment 
for the production of new cloud line. 

C. THESIS OBJECTIVE 

The objective of this study is to observe and describe the 
spatial and radiative variations of ship tracks that form in 
the clouds associated with transition regions. Chapter II will 
outline the approach of the study and describe the data 
collection and analysis. Chapter III will discuss the results 
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and Chapter IV will present the conclusions and make 
recommendations for future study. 
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II. APPROACH 



A. OVERVIEW 

The studies noted earlier clearly indicate variations in 
the liquid water content, vertical velocity, and inversion 
height in the ambient clouds associated with cloudiness 
transition regions (CTR) . The focus of this study is to locate 
ship tracks that transit these regions and describe their 
radiative and spatial characteristics as a function of the 
surrounding regime. For the purpose of this study the ambient 
clouds associated with CTR's are divided into cumulus, broken 
stratocumulus , and solid stratus similar to the method 
utilized by Betts and Boers, 1990. Since in situ measurements 
of liquid water content, vertical velocity and inversion 
height are not available for study, mean quantities for each 
cloud type will be based on the data collected in previous 
studies of marine stratiform clouds by Brost et al, 1982. 
These values will be utilized to characterize the CTR 
environments and provide some insight into the physical 
mechanisms responsible for ship track formation. 
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B . DATA 



The data utilized in this study are from the Advanced Very 
High Resolution Radiometer (AVHRR) onboard the NOAA-9/10 
satellites. These satellites are in a polar orbit 
approximately 850 kilometers above the earth's surface. The 
instrument continuously records 2048 samples per scan line 
centered on the nadir. It scans upwelling radiation, both 
emitted and reflected energy, over five wavelength channels 
centered at 0.63, 0.86, 3.7, 11, and 12 ^m. This scanning 
geometry produces a pixel resolution of 1 km by 1 km at nadir. 

The primary data source utilized in this study come from 
the First ISSP Regional Experiment (FIRE) which was conducted 
from 28 June to 14 July 1987 off the coast of southern 
California. Both morning and afternoon passes were collected, 
copies of which are archived in the Naval Postgraduate School 
IDEA Laboratory. Using the IDEA Lab VAX computer system, the 
AVIAN software package was employed to glean and navigate the 
magnetic tapes ^litaining these data. This allowed the 
production of an overview of the entire satellite pass. This 
overview is visually scanned for the presence of ship tracks. 
When ship tracks are found in CTR's or cloud regimes 
associated with them, a 512 km x 512 km channel 3 subscene is 
extracted from the overview for further analysis and ship 
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track digitization. After the track is digitized, five 
additional image products are also extracted from the 
overview. A brief description of the characteristics of these 
image products is provided below: 



LOWl - 



L0W3 - 



S12A - 



TMP4 - 
T45 - 



channel 1 albedo scaled by the solar zenith 
angle and low cloud asymmetric reflectance 
factor (in percent) . 

channel 3 albedo scaled by the solar zenith 
angle and low cloud anisotropic reflectance 
factor (in percent) . 

a ratio of channel 1 albedo to channel 2 
albedo (1.0 to 3.0). 

channel 4 brightness temperature (°K) 
difference between channel 4 and channel 5 
brightness temperatures (°K) . 



These 512 X 512 image products subscenes become input 
parameters into the ship track data extraction algorithm to te 
described next. 
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C. SHIP TRACK DATA EXTRACTION ALGORITHM 



Following on the approach of Coakley et al . (1987) and 
Morehead (1988) , a ship track detection algorithm was 
developed at the Naval Postgraduate School. A description of 
the algorithm is contained in Nielsen and Durkee, (1992) . The 
algorithm divides the subscene into 16 Ian x 16 km subareas and 
computes the mean and standard deviation for each subarea. 
Next, it performs a brightness "neighborhood" test to 
determine which local pixels show some relation to other local 
pixels. Then the algorithm determines which brightness 
neighborhoods belong to ship tracks and constructs a linear 
pathway connecting the neighborhoods. 

The algorithm was designed to automate the ship track 
detection process by independently identifying ship tracks 
with a minimum of manual intervention. In this study the 
algorithm was modified to allow the user to manually select 
the ship track (s) of interest, digitize a series of points to 
define the track and input the points into the algorithm. This 
greatly speeds up the extraction process since the algorithm 
only scans the subscene for the specified track and ignores 
other tracks and anomalous linear features. 
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Once the track is selected, digitized, and inputted into 
the algorithm, it extracts a 31 km swath of data along the 
entire length of the track. This swath of data contain the 
radiative signature of both the ship track and the ambient 
cloud field surrounding it. A file is created for each of the 
subscenes previously described. It is these data which will be 
analyzed to determine the variation in the spatial and 
radiative characteristics of ship tracks in CTR's. 

D. DATA ANALYSIS 

Before any averaging technique or statistical analysis can 
be performed, the data contained in the linerized ship track 
records file must be transferred to 1 km x 1 km grid to ensure 
that the size of a pixel is constant throughout the width and 
length of the track. This adjustment is necessary because the 
combined effect of the AVHRR's scan geometry and the earth's 
curvature causes a degradation in resolution. Additionally, 
the distance between data points increases away from the nadir 
of the satellite, particularly near the edges. So the 
extracted swath of data contain some holes particularly when 
curved tracks are straightened out. The best resolution is 
obtained directly below the satellite, at nadir, resulting in 
a pixel size of approximately 1.1 km x 1.1 km spaced 
approximately .9 km apart. This resolution decreases toward 
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either end of the scan line where the pixels measure about 2.4 
km along- track and about 14.7 km across- track. Basically, the 
algorithm takes the linerized records file, determines its 
location in the overview, and applies a scaling correction to 
each pixel based on the angular distortion present at that 
point. Then the records file is mapped onto a 1 km grid. Each 
data point of the records file is assigned an x-y coordinate 
in the grid. Values are assigned to the grid points based on 
a weighting function which assigns a 1/r^ weight to each data 
point within a user defined radius of 2-5 km around each grid 
point. This variable radius is necessary because the distance 
between data points varies as a function of the tracks' 
location on the overview. It also allows the user some control 
over the number of data points used in computing each grid 
point value. The procedure results in some smoothing but the 
ship track/ambient cloud swath essentially retains its 
radiative signature and spatial characteristics. A raw 
satellite image and a linerized, weighted-average version of 
the same ship track are shown in Figure 4 . 

E . AVERAGING METHODS 

The procedure employed in the previous section allows 
ship tracks to be averaged and statistically analyzed. Sample 
pixels are extracted from the records file for each subscene. 
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Figure 4. Raw and Linearized, Weighted Average Versions of 
the Same Ship Track. 
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The extraction algorithm takes the pixel at the center ship 
track and one pixel on either side. Then, using a brightness 
gradient to find the edge of the ship track, it extracts six 
pixels from the uncontaminated cloud, three from either side 
of the track. A distance of 15 km from the edge is utilized as 
a separation distance. The track pixels are averaged together 
and the six cloud pixels are averaged together for each line 
of data along the track. This operation is performed for each 
line of data in the records file. Figure 5, which is an 
enlargement of the previous weight -averaged ship track, 
depicts pixel scale data points utilized for the ambient cloud 
and ship track data sets. Combining records files for similar 
track/cloud types allows class averages and comparisons to be 
produced for each ship track/cloud pair. 

Once the tracks for each cloud type are averaged, the 
spatial and radiative variations are examined. A summary of 
the statistics utilized in the study for each ship track/cloud 
pair are summarized in Table I. 
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Figure 5. Typical Ship Track and Ambient Cloud Pixels 
Utilized in Averaging Scheme, 



The results of the analysis of the radiative and spatial 
variations of ship tracks, with respect to the environment in 
which the form, are detailed in Chapter III. 
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Table I SUMMARY OF STATISTICAL DATA COLLECTED ON SHIP 
TRACKS AND AMBIENT CLOUD FIELDS. 



LOWl AVERAGE AMBIENT CLOUD ALBEDO 

LOWl AVERAGE TRACK ALBEDO 

LOW3 AVERAGE AMBIENT CLOUD ALBEDO 

LOW3 AVERAGE TRACK ALBEDO 

S12A AVERAGE AMBIENT CLOUD ALBEDO 

S12A AVERAGE TRACK ALBEDO 

TMP4 AVERAGE AMBIENT CLOUD TEMPERATURE 

TMP4 AVERAGE TRACK TEMPERATURE 

T45 AVERAGE AMBIENT CLOUD TEMPERATURE DIFFERENCE 
T45 AVERAGE TRACK TEMPERATURE DIFFERENCE 
DELTA PERCENT CHANGE LOWl CLOUD/TRACK 
DELTA PERCENT CHANGE LOW3 CLOUD/TRACK 
DELTA PERCENT CHANGE S12A CLOUD/TRACK 
DELTA PERCENT CHANGE TMP4 CLOUD/TRACK 
DELTA PERCENT CHANGE T45 CLOUD/TRACK 
WIDTH OF TRACK VS. LENGTH 
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III. RESULTS 



A. OVERVIEW 

A total of 63 ship tracks were collected and analyzed in 
this study. They were collected from both morning and 
afternoon satellite passes. They had a latitudinal range from 
approximately 20 °N to 45 °N and a longitudinal range from 118 °W 
to 138°W. Throughout the period of the study the entire region 
was under the influence of the semi -permanent subtropical 
high. Surface pressures ranged from 1012 mb to 1024 mb. 
Surface winds were generally north-northwesterly varying from 
5 to 30 kt throughout the period. 

Ship tracks were observed in a wide range of cloud types, 
textures and patterns. These were analyzed and classified into 
three broad categories: solid stratiform, broken stratiform 
and cxomulus . The observed radiative and spatial variation by 
satellite image product and cloud type are described by 
analyzing the class averages for the statistical data 
described previously. The results of this analysis are 
summarized in the following sections. 
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B. STATISTICAL METHODS AND TESTS 



The analysis of these data required the determination of 
the mean and standard deviation for the track and ambient 
cloud for each cloud type and image product. The means of the 
track and cloud values were compared and a percentage assigned 
based on the differences observed between the means. An 
additional statistic was computed by subtracting the value of 
each ambient cloud pixel from the value of its corresponding 
track pixel. This difference was divided by the value of the 
ambient cloud pixel to produce a delta percent change (DPC) 
between the track and cloud. The mean and standard deviation 
of average of the differences were also computed. 

The Mann-Whitney (MW) ranking scheme for two independent 
samples, Snedecor and Cochran (1967) , was utilized to test the 
significance of the difference between the means of the track 
and cloud values. The null hypothesis that both samples come 
from the same population is rejected if the MW statistic gives 
a probability value that produces a significance level greater 
than 95 percent. In each of the following cases the MW test 
results will be reported to estimate the probability that the 
the differences observed between two means is significant and 
not the result of chance. 
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C. DATA DISPLAY METHOD 



Prior to any discussion of the results of this study, a 
description of data display methods is necessary. The data 
collected for each cloud/track pair within a particular cloud 
regime are displayed in both scatterplot and histogram form. 
Both displays depict the cloud regime or track with the head 
of the track to the left of the graph. Figure 6 illustrates 
scatterplots for a single ship track in both LOWl and L0W3 . 
Combining a number of tracks produces the plots contained in 
the following section. Although there is considerable 
variability in one ship track, most of the variation present 
is from track to track. 

D. OBSERVATIONS OF RADIATIVE CHARACTERISTICS 
1. Stratus Clouds and Associated Tracks 
a. LOWl 

The LOWl mean for the albedo of stratus clouds is 
34.672 with a standard deviation of 13.919 while the mean for 
the tracks which form in them is 38.652 with a standard 
deviation of 14.512. This indicates that, on the average, 
tracks that form the stratus cloud are 11% brighter than the 
ambient cloud. Utilizing the MW test, the probability that the 
means are from different populations is 100%. The mean of the 
DPC between track and cloud is 15.128 with a standard 
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LOWl ST TRACK ALBEDO 




LOW3 ST TRACK ALBEDO 




Figure 6. LOWl and LOW3 Scatterplots Depicting a Single Ship 
Track. 
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deviation of 31.004. Figure 7 shows the raw data utilized to 
compute the means. The LOWl albedo value for stratus range 
from 8 to 77 with 23.1 being the most frequently occurring 
value. The cloud brightness for individual track is fairly 
constant along each track but varies considerably between 
tracks. Figure 8 shows the albedo values for ship tracks in 
stratus. The values vary from 10 to 78 with a value of 23.8 
occurring most frequently. The delta percent change (DPC) 
shown in Figure 9 fluctuates between -30% and 60% with a 
maximum occurring value of 6%. These means values are 
consistent with, but higher than, the observations reported by 
Coakley et al . , 1987, who noted that the change was caused by 
an increase in cloud droplets numbers. A larger concentration 
of cloud droplets causes the cloud to be more reflective, 
jb. LOW3 

The L0W3 mean for the albedo of stratus clouds is 
8.191 with a standard deviation of 5.105. The mean for the 
ship tracks in them is 10.397 with a standard deviation of 
5.643. The percent change in reflectivity on the average 
between the track and the cloud is 27%. The MW shows a 100% 
probability the difference is significant. The mean of the DPC 
is 35.631 with a standard deviation of 40.55. Figure 10 
displays the plots for L0W3 stratus clouds. The values vary 
from 2 to 25 with maximum occurring near 4.1. The ship track 
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LOWl AVG ST CLOUD ALBEDO 
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Figure 8. Scatterplot and Histogram for LOWl Stratus Tracks. 
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data points are displayed in Figure 11. They range from 4 to 
30 with a most frequent value at 5.1. The DPC, shown in 
Figure 12, ranges between -30% and 260% but the most 
frequently occurring value is 16.8%. The injection of 
additional CCN into the cloud by a passing ship causes smaller 
water droplets to form. Since the smaller in- track water 
droplets are more efficient scatters of energy than the 
neighboring larger ambient cloud droplets, ship tracks tend to 
be more reflective at 3.7 /im (Coakley et al., 1987). 

C. S12A 

The mean of the S12A ratio for stratus clouds is 
1.175 with a standard deviation of 0.008. The ship track class 
average is 1.174 with a standard deviation of 0.008. This 
states that both cloud and track are brighter in channel 1 
than channel 2 . On the average the percent change between 
cloud and track is slightly negative at -0.004. The MW shows 
a 79% probability the difference is significant. The mean of 
the DPC is -0.003 with a standard deviation of 2.210. These 
observations are consistent with the different channel 1 and 
2 wavelengths. Both the ambient cloud and the ship track are 
brighter in Channel 2. However, because of the shift to 
smaller, more reflective water droplets in the track, ship 
tracks gets brighter faster than the ambient clouds. This same 
effect is observed at 3.7 ^m where ship tracks are 
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Regime . 
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significantly brighter than the ambient cloud. The resultant 
ratio of the S12A for the ship track and cloud is therefore 
negative. The scatterplots and histograms for the data show 
little variation. For completeness, the S12A plots for 
stratus, broken stratus and cumulus are included in the 
Appendix. 

d. TMP4 

The mean for TMP4 temperature in both stratus 
cloud and the ship track is 285.2 °K. The standard deviations 
were 2.358 °K and 2.412 °K, respectively. The MW shows a 56% 
probability the difference is significant. The mean of the DPC 
is -0.001 with a standard deviation of 0.173. These means are 
identical indicating that both the cloud and track emit 
equally as well at 11 /xm. It also suggests that ships cause no 
appreciable change in the cloud top height or a temperature 
difference would be observed. The scatterplots and histograms 
for the data show little variation between cloud and track. 
The TMP4 plots for stratus, broken stratus and cumulus are 
also presented in the Appendix. 

e. T45 

The mean for the T45 temperature difference is 
-.395 °K with and standard deviation of 1.306 °K while the 
mean for ship tracks is 0.001 °K with a standard deviation of 
.306 °K. This translates to a large percent change from track 
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to cloud but represents the differences between numbers on the 
order IE- 02 and IE- 03. The sensor can only resolve temperature 
differences on the order of IE- 01. So, although a weak 
signature appears to exist no real significance is attached to 
the large percent change observed. The scatterplots and 
histograms for this data, broken stratus and cumulus are 
provided in the Appendix. 

A summary of the statistics change for each of the image 
products described above for stratus clouds is contained in 
Table II. 

2. Broken Clouds and Associated Tracks 
a. LOWl 

The LOWl mean albedo of broken stratus clouds and 
ship tracks is 23.422 and 27.770 respectively. The standard 
deviation for the cloud average is 8.535 and versus 9.039 for 
the track. On the average, the total increase in albedo 
between the track and the cloud is 18%. The MW test shows this 
to be significant to 100%. The mean of the DPC is 26.016 with 
a standard deviation of 44.854. The average percent change and 
the DPC are both higher than that observed in the stratus 
case. In broken clouds the ship track is contrasted against a 
background clouds that have gaps. The radiative signature of 
these clouds is influenced the less reflective sea surface. 
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Table II. SUMVIARY OF STATISTICS FOR STRATUS CLOUDS /TRACKS . 



PRODUCT 


MEAN 


STD DEV 


% CHANGE 


LOWl ST 


34.652 (A) 


13.919 




LOWl TRACK 


38.341 (A) 


14.512 




LOWl TRK/ST 






0.11 


LOWl DPC 


15.128 (A) 


31.004 




LOW3 ST 


8.191 (A) 


5.105 




LOW3 TRACK 


10.397 (A) 


5.643 




LOW3 TRK/ST 






0.27 


LOW3 DPC 


35.631 (A) 


40.55 




S12A ST 


1.176 (R) 


0.008 




S12A TRACK 


1.174 (R) 


0.008 




S12A TRK/ST 






-0.001 


S12A DPC 


-0.003 (R) 


2.210 




TMP4 ST 


282.2 (T) 


2.358 




TMP4 TRACK 


282.2 (T) 


2.412 




TMP4 TRK/ST 






0.000 


TMP4 DPC 


-0.001 (T) 


0.173 




T45 ST 


-0.396 (D) 


1.306 




T45 TRACK 


0.000 (D) 


0.f08 




T45 TRK/ST 






-100.00 


T45 DPC 


-9.280 (D) 


195.47 




A=ALBEDO, R=RATIO, T=TEMPERATURE (®K), D=DIFFERENCE 
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This makes ship tracks which form in broken clouds appear 
brighter relative to their environment. The albedo values for 
LOW 1 broken stratus vary from 7 to 58 with a most frequently 
occurring value of 20 as shown in Figure 13. For the track. 
Figure 14 shows the values fluctuate between 9 and 64 the 
maximum occurring near 24. Figure 15 illustrates the DPC 
between track and cloud. It ranges from -60% to 237% with a 
maximum at 8%. As noted the tracks in broken cloud have a 
lower albedo than in stratus. One possible explanation for 
this is that the ambient broken clouds have a lower liquid 
water content (LWC) . The additional flux of CCN into the cloud 
causes a shift to smaller droplet size but because the LWC in 
broken clouds is less than that in stratus clouds, the total 
concentration of smaller radii water droplets is lower. 
b. LOW3 

The mean for LOW3 broken cloud albedo is 5.842 
with a standard deviation of 3.114. The track mean albedo is 
8.210 with a standard deviation of 3.820. On the average there 
is a 40% increase in reflectivity between the track and the 
cloud. The MW test shows the significance of the differences 
in the means to be significant to 100%. The mean of the DPC is 
53.911 with a standard deviation of 57.644. Figure 16 presents 
the LOW3 broken stratus cloud data. The values range from 8 to 
15 with bimodal maximums occurring at 2.8 and 7.7. This 
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Regime . 
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bimodal signature results from the variability observed 
between tracks. The track values, Figure 17, fluctuate between 
1.6 and 18 with a most frequently occurring value at 8.6. The 
DPC, shown in Figure 18 has a wide range with a minimum near 
14.3%. While this is substantial increase over the stratus 
case, it is consistent with smaller droplet radii and 
increased droplet concentration over the ambient cloud. It is 
also consistent with the comparable increase observed in LOWl. 



C. S12A 

The S12A mean ratio is 1.145 with a standard 
deviation of .102 while the mean for the ship track is 1.136 
with a standard deviation of 0.009. The average percent change 
between the track and the cloud is -0.008. The MW test places 
the significance of the difference at 100%. The mean of the 
DPC is -0.669 with a standard deviation of 3.285. The same 
effect previously discussed in stratus is present in broken 
clouds as well, albeit more pronounced. Since the contrast 
between broken clouds and ship tracks is greater the resultant 
S12A ratio of track to cloud is slightly more negative than 
that observed in stratus. Scatterplots and histograms from 
which these data are derived are in the Appendix. 
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Figure 16. Scatterplot and Histogram for L0W3 Broken Stratus 
Regime . 
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Figure 17 . Scatterplot and Histogram for L0W3 Broken Stratus 
Tracks . 
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d. TMP4 



The mean for TMP4 temperature in broken cloud is 
284.4 °K with a standard deviation of 1.921 °K. The mean for 
the track is 284.3 °K with a standard deviation of 1.881 °K. 
The MW significance for the difference in the means is 99%, 
much higher that in the stratus case. The mean of the DPC is 
-0.005 with a standard deviation of 0.025. The influence of 
the sea surface temperature is perhaps one reason these 
statistics are different the stratus case. The statistics 
point to a difference in the means but is detectable 
difference exists in the ambient cloud and track temperature. 
The Appendix contains the scatterplots and histograms for the 
above data . 

e. T45 

The mean T45 temperature difference is -0.261 °K 
for broken stratus with a standard deviation of 1.417 °K. The 
average difference for the track is 0.200 °K with a standard 
deviation of 0.343 °K. The percent change is again large but 
not considered significant. Scatterplots and histograms for 
this data are included in the Appendix. 

A summary of the statistics each of the image products 
described above for broken stratus clouds is contained in 
Table III. 
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Table III. SUMMARY OF STATISTICS FOR BROKEN STRATUS 
CLOUD /TRACKS . 



1 PRODUCT 


MEAN 


STD DEV 


% CHANGE 


LOWl BK 


2S.422 (A) 


8 .5S5 




LOWl TRACK 


27.770 (A) 


9 .0S9 




LOWl TRK/BK 






0.18 


LOWl DPC 


26 . 016 


44.854 




LOWS BK 


5.842 (A) 


2.114 




LOWS TRACK 


8.210 (A) 


S .820 




LOWS TRK/BK 






0.40 


LOWS DPC 


5S .911 


57.644 




S12A BK 


1.145 (R) 


0.102 




S12A TRACK 


1.1S6 (R) 


9 .0S9 




S12A TRK/BK 






-0.008 


S12A DPC 


-0.669 


S.285 




TMP4 BK 


284.4 (T) 


1.921 




TMP4 TRACK 


284.2 (T) 


1.881 




TMP4 TRK/BK 






0.000 


TMP4 DPC 


-0.005 


B.02# 




T45 BK 


-0.261 (D) 


1.417 




T45 TRACK 


0.200 (D) 


0.242 




T45 TRK/BK 






-176.00 


T45 DPC 


-20.97 


182.81 




A=ALBEDO, R=RATIO, T=TEMPERATURE (°K), D=DIFFERENCE 
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3. Ciimulus Clouds and Associated Tracks 



a. LOWl 

The mean of the albedo of cumulus clouds is 
19.054 with a standard deviation of 8.210. The mean albedo of 
the ship tracks that form in them is 23.42 with a standard 
deviation of 9.464. On the average there is a 22% increase in 
albedo between the track and cloud which implies that cumulus 
tracks are much brighter than both stratus tracks and broken 
stratus tracks relative to their environment. The Mann- 
Whitney test gives a 100% probability the difference between 
these two means is significance. The mean of the DPC is 
30.231 with a standard deviation of 46.975. Figure 19 shows 
the data distribution within the LOWl cumulus cloud regime. 
There is a range of values from 5 to 47 with a maximum 
occuring at 13.5. The track albedo data points, shown in 
Figure 20, vary from 7 to 52 with the most frequently 
occurring value at 13.4. Figure 21 show the DPC for LOWl 
between track and cloud. The values fluctuate over a wide 
range with a maximum occurring a 10.2%. The radiative 
signature of cumulus clouds is heavily influenced by the low 
reflectivity of the sea surface. This cloud type also has the 
lowest LWC of all three cases because the inversion is weaker 
and allows the cumulus clouds to penetrate and disperse 
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moisture. In this type of cloud the ship usually produces new 
cloud as it forms a track. 
b. LOW3 

The mean for the L0W3 albedo in cumulus clouds is 
5.748 with a standard deviation of 2.945. The mean for the 
tracks is 7.721 with a standard deviation of 2.864. The 
average percent change in albedo is 34%. The MW test assigns 
a significance of 100% to the difference between the means. 
The data for L0W3 ambient cumulus cloud regime is shown in 
Figure 22. They vary from 1.2 to 12 with a maximum value at 
4.5. Figure 22, which displays the track data, has a range 
from 2 to 15 with a maximum of 6.6. The DPC for L0W3 cumulus 
track and cloud shows a large fluctuation with a maximum at 
13.7%. This is a substantial increase in albedo over the 
ambient cloud but reverses an increasing trend from stratus to 
broken. This is consistent with a lower LWC in cumulus clouds 
due to vertical dispersion. That implies that there is less 
total water droplet to be affected by the injection of CCN so 
while a shift to smaller droplet size occurs it is not a 
dramatic as the shift in stratus and broken stratus clouds. 
The mean of the DPC is 51.778 with a standard deviation of 
61.425. 
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Figure 19 . Scatterplot and Histogram for LOWl Cumulus Cloud 
Regime . 
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C. S12A 



The S12A mean ratio is 1.172 in cumulus clouds 
with a standard deviation of 0.007 while the mean ratio of the 
track is 1.158 again with a standard deviation of 0.007. The 
average total percent change is -0.011 is larger than the 
previous two cases and establishes a trend that ship tracks 
gets brighter relative to their environment as the cloud 
fraction decreases. The MW test gives a significance of 100% 
to the differences in the means. The mean of the DPC is -1.088 
with a standard deviation of 2.676. Scatterplots and 
histograms for the above data are contained in the Appendix. 

d. TMP4. 

The mean for the TMP4 for cumulus clouds is 286.5 
°K and 286.4 °K in the track. There is no appreciable 
difference between the means of the cloud and track 
temperatures. However, the MW test gives a significance to the 
difference of 100%. The mean of the DPC is -0.005 with a 
standard deviation of 0.136. The cloud and track are both 
warmer than both the stratus and broken cases which is 
consistent with lower clouds. The Appendix includes the 
scatterplots and histograms for this data. 

e. T45 

The mean of the T45 difference for cumulus clouds 
is -0.154 °K with a standard deviation of 1.236 °K. The track 



51 




52 




53 




Track/Cloud. 
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mean is 0.228 °K with a standard deviation of 0.233 °K. This 
again represents a large percentage change but no statistical 
significance is attached to this difference. Scatterplots and 
histograms for this data are in the Appendix. 

A summary of the statistics for each of the image products 
described above for cumulus clouds is contained in Table IV. 
Table V contains a summary of the Mann -Whitney test results. 

E. OBSERVATIONS OF SPATIAL CHARACTERISTICS 

Each of the 63 ship tracks collected for the study was 
manually analyzed to determine the width of track as a 
function of length. The tracks ranged in length from tens to 
hundreds of kilometers. One ship track observed, but not 
analyzed in this study, was well over 1000 km long. If the 
ship was moving at 25 knots then this track was over two days 
old. 

1. Width of Tracks 

The spatial variations between cloud types was not as 
dramatic as originally hypothesized. One hundred kilometers 
down the length of the track was used as a benchmark for the 
analysis of track width and dispersion. Tracks which form in 
stratus cloud were found to be 4-9 km wide. Broken stratus 
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Table IV. SUMMARY OF STATISTICS CUMULUS CLOUD/TRACK. 



PRODUCT 


MEAN 


STD DEV 


% CHANGE 


LOWl CU 


19.054 (A) 


8.210 




LOWl TRACK 


23.427 (A) 


9.564 




LOWl TRK/CU 






0.22 


LOWl DPC 


30.231 


46.975 




LOW3 CU 


5.748 (A) 


2.945 




LOW3 TRACK 


7.721 (A) 


2.864 




LOW3 TRK/CU 






0.34 


LOW3 DPC 


51.778 


61.425 




S12A CU 


1.172 (R) 


0.007 




S12A TRACK 


1.158 (R) 


0.007 




S12A TRK/CU 






-0.011 


S12A DPC 


-1.088 


2.676 




TMP4 CU 


286.5 (T) 


1.443 




TMP4 TRACK 


286.4 (T) 


1.431 




TMP4 TRK/CU 






.000 


TMP4 DPC 


-0.005 


0.136 




T45 CU 


-0.154 (D) 


1.236 




T45 TRACK 


0.228 (D) 


0.233 




T45 TRK/CU 






-248.00 


T45 DPC 


-17.233 


120.286 




A=ALBEDO, R=RATIO, T=TEMPERATURE (°K), D=DIFFERENCE 
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Table V. SUMMARY OF MANN-WHITNEY TEST RESULTS. 



MEAN VALUES CLODD/TRACK 


PROBABILITY 


LOWl STRATUS /TRACK 


100% 


LOWl BROKEN/TRACK 


100% 


LOWl CUMULUS /TRACK 


100% 


LOW3 STRATUS /TRACK 


100% 


LOW3 BROKEN/TRACK 


100% 


LOW3 CUMULUS /TRACK 


100% 


S12A STRATUS /TRACK 


>79% 


S12A BROKEN/TRACK 


100% 


S12A CUMULUS /TRACK 


100% 


TMP4 STRATUS /TRACK 


>56% 


TMP4 BROKEN/TRACK 


>99% 


TMP4 CUMULUS /TRACK 


>99% 


T45 STRATUS /TRACK 


100% 


T45 BROKEN/ TRACK 


>93% 


T45 CUMULUS /TRACK 


100% 



cloud regime tracks ranged from 4-10 km. The width of tracks 
in cumulus clouds at 100 km also was between 4-9 km. This 
result was surprising since the weaker inversion, particularly 
in cumulus clouds, allows vertical motion to penetrate and 
disperse boundary constituents. It was expected ship tracks in 
the cumulus environment would be wider and more diffuse than 
was observed. 
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2 . Length of Tracks 

The length of ship tracks in each environment was 
difficult to quantify. Tracks were collected in each image 
product so as to include the head of the track and as much of 
the track as would fit in the 512 X 512 image. Some tracks 
were too long to be captured in one image. In addition, the 
ship track extraction algorithm limited extraction of tracks 
too near to the edges of the subscene to ensure a complete 
sample of the track and the ambient cloud. So any quantitative 
analysis of ship track length in this study would not be 
valid. Qualitatively however, observations throughout the 
collection and analysis process indicate that tracks which 
form in stratus clouds tend to be longer. This is consistent 
with the strong capping inversion associated with this 
environment. This cap would tend to limit vertical dispersion 
of the track and perhaps make it more long lived. Broken and 
cumulus tracks were generally the same, although exceptionally 
long tracks were observed in these environments as well. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The purpose of this thesis was to investigate and analyze 
the radiative and spatial variations of ship tracks in three 
distinct marine stratiform cloud environments: stratus, broken 
stratus and cumulus. The variations is ship tracks are 
consistent with the differences in inversion strength, liquid 
water content and vertical velocity present in each 
environment . 



The results indicate that the radiative signature of ship 
tracks varies with respect to the environment in which it 
forms. Specifically; 



• A comparison of the means of ship tracks and ambient 
clouds for each environment shows that the tracks are at 
least 11% brighter in LOWl and at least 27% brighter in 
L0W3 . Cumulus tracks have the highest average albedo 
relative to their environment in LOWl at 22%. Broken 
stratus tracks are highest in their environment at 40% 
in LOW3. 

• The delta percent change statistic indicates broken 
stratus tracks are 2% brighter than cumulus tracks 
and 18% brighter than stratus tracks relative to 
their environment in L0W3 . In LOWl cumulus tracks are 
4% brighter than broken stratus tracks and 15% 
brighter than stratus tracks. 
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There is a link between the albedo of the track 
and the liquid water content of the cloud. Ctunulus 
tracks are generally brightest with respect to their 
environment but lowest in absolute albedo over stratus 
and broken stratus tracks. Ship do not appear to enhance 
the liquid water content of the clouds they pass 
beneath. 

The S12A signature observed supports the observations 
noted in LOWl and L0W3 . 

There was no temperature change observed between the 
track and its environment in all three cases. 

The T45 signature shows a large percentage change 
between the track and the ambient cloud but the 
differences are between very small numbers, so the 
percent change is greatly inflated. 

Spatial variations in the width of ship tracks between 
each environments was small but the lengths of tracks in 
solid stratus tended to be longer than in the other two 
environments . 



B . RECOMMENDATIONS 

This study represents the first attempt at identifying 
differences in the radiative and spatial characteristics of 
ship tracks with respect to the cloud type in which they form. 
In this approach cloud fields were classified into three broad 
categories. Many cloud sub -types were observed within each 
broad category. The physical mechanisms responsible for ship 
track formation are still not well understood so future ship 
track research should: 

• Continue to refine the ship track detection algorithm. 
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• Automate the collection, extraction, classification and 
archiving of ship tracks . 

• Continue to analyze ship tracks in a variety of 
cloud environments for clues to their formation 
mechanisms . 

• Continue in-situ location of ships producing tracks 

• Continue to investigate ship track formation mechanisms 
with respect to cloud liquid water content and CCN 
concentrations to better understand the environmental 
parameters necessary for ship track production. 



Much work remains to be done on this subject. Collection 
and analysis of this data is important because of its 
applications to military intelligence, commerce, and 
climatology. Continuing investigation of these small cloud 
streaks could, perhaps, help unravel the effect man has on the 
global environment. 
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APPENDIX. SCATTERPLOTS AND HISTOGRAMS 



The scatterplots and histograms for S12A, TMP4 and T45 
stratus, broken stratus and cumulus cloud regimes and ship 
tracks can be found in this Appendix. The S12A plots appear in 
Figures A1 to A9 . The TMP4 plots appear in Figures AlO to A18 
and the T45 plots appear in Figures A19 to A27. 
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Figure Al. Scatterplot and Histogram for S12A Stratus Cloud 
Regime . 
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Figure A2 . Scatterplot and Histogram for S12A Stratus Track. 



64 



w 

o 

<c 

£- 

Z 

H 

CJ 

crS 

Cu 



S12A PERCENT CHANGE TRACK/ST CLOUDS 

10 



0 



TRACK LENGTH (KM) 

PeRCENT OCCURRENCE 

t.^ 13.5 20^2 ^.0 

10.8 



13. <5 



11.0 



<9.2<? 

LTi 

^ 0 <J8 

I ..» 



LTi 

</l 

3.333E-01 

-2.28 



4 60 



-6 <J1 



-9.22 
-11 .5 

O00 1200 1800 2400 

FREQUENC'^ 





Figure A3. Scatterplot and Histogram for DPC S12A Stratus 
Track/Cloud* 
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Figure A4. Scatterplot and Histogram for S12A Broken Stratus 
Cloud Regime. 
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Figure A6. Scatterplot and Histogram for DPC S12A Broken 
Stratus Track/Cloud. 
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Figure A7 . Scatterplot and Histogram for S12A Cumulus Cloud 
Regime . 
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Figure A8. Scatterplot and Histogram for S12A 
Tracks . 
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Scatterplot and Histogram for TMP4 Stratus Cloud 
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Figure All. Scatterplot and Histogram for TMP4 Stratus 
Tracks . 
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Figure A12 . Scatterplot and Histogram for DPC TMP4 
Track/Cloud. 
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Figure A13. Scatterplot and Histogram for TMP4 Broken 
Stratus Clouds. 
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Figure A14 . Scatterplot and Histogram for TMP4 Broken 
Stratus Tracks. 
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Figure A15. Scatterplot and Histogram DPC TMP4 Broken 
Stratus Track/Cloud. 
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Figure A16. ScatCerploC and Histogram for TMP4 Cumulus Cloud 
Regime . 
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Figure A17 . Scatterplot and Histogram for TMP4 Cumulus 
Tracks . 
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Figure A18. Scatterplot and Histogram for DPC TMP4 Cumulus 
Track/Cloud. 
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Figure A19 • Scatterplot and Histogram for T45 Stratus Cloud 
Regime . 
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Figure A21. Scatterplot and Histogram for DPC T45 Stratus 
Track/Cloud, 
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Figure A22. Scatterplot and Histogram for T45 Broken Stratus 
Cloud Regime. 
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Figure A23 . Scatterplot and Histogram for T45 Broken Stratus 
Tracks . 
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Figure A24. Scatterplot and Histogram for DPC T45 Broken 
Stratus Track/Cloud. 
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Scatterplot and Histogram for T45 Cvimulus Cloud 
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Figure A26. Scatterplot and Histogram for T45 Cumulus 
Tracks . 
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Figure A27. Scatterplot and Histogram for DPC T45 Cumulus 
Track/Cloud. 
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